Phonon Thermodynamics versus Electron-Phonon Models 
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Applying the path integral formalism we study the equilibrium thermodynamics of the phonon 
field both in the Holstein and in the Su-Schrieffer-Heeger models. The anharmonic cumulant series, 
^ ' dependent on the peculiar source currents of the e-ph models, have been computed versus tempera- 

ture in the case of a low energy oscillator. The cutoff in the series expansion has been determined, 
in the low T limit, using the constraint of the third law of thermodynamics. In the Holstein model, 
the free energy derivatives do not show any contribution ascribable to e-ph anharmonic effect. We 
find signatures of large e-ph anharmonicities in the Su-Schrieffer-Heeger model mainly visible in the 
temperature dependent peak displayed by the phonon heat capacity. 
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I. INTRODUCTION 
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^ ■ The equilibrium thermodynamical properties of real 
systems provide signatures of phonon anharmonicities in- 
3 ■ duced both by electron-phonon and phonon-phonon in- 
teractions. While in general the exact thermodynamics 
of anharmonic e-ph model Hamiltonians can be derived 
only in one dimension and in the classical limit, vari- 
ational methods based on the Bogoliubov's inequality 
• [1] permit to obtain i) upper bounds for the partition 
functions [2] and ii) a correct estimate of the entropy 
within the harmonic theory once the harmonic phonons 
are self-consistently replaced by the temperature depen- 
. dent phonon frequencies [3] . Such a replacement however 
' does not hold for other thermodynamic functions and 
constant volume anharmonic effects have to be explicitly 
computed by diagrammatic perturbative methods [4,5]. 
As at high temperatures phonon-phonon anharmonici- 
ties are large, in the intermediate to low temperature 
range the e-ph coupling strength may induce the dom- 
inant anharmonic effects in the phonon subsystem, the 
latter being dependent on the peculiarities of the e-ph 
model Hamiltonian. 

The Holstein Hamiltonian [6] and the Su-Schriffcr- 
Heeger Hamiltonian [7] have been widely investigated in 
the last decades also in view of the growing interest for 
polaronic mechanisms in high T c superconductors [8] and 
polymers [9]. While the Holstein model assumes a local 
coupling of tight binding electrons to optical phonons, 
in the SSH model the electron-phonon interaction mod- 
ifies the electron hopping matrix elements thus leading 
to a non-local (in momentum space) coupling with ver- 
tex function dependent on both the electron and the 
phonon wave vector [10,11]. As the ground state of the 
SSH Hamiltonian is twofold degenerate localized solitonic 
solutions appear in the system together with polarons 
[12,13] whose formation depends however on the strength 
of the electron coupling to the induced lattice deforma- 
tion [14] and on the value of the adiabatic parameter [15]. 
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Analysis of the Holstein and SSH Hamiltonians gen- 
erally assume harmonic phonon models although lattice 
anharmonicities may induce relevant microscopic effects 
such as the reduction of the Holstein polaron effective 
mass [f 6] along with a consistent broadening of the po- 
laron size [17]. In large polaron transport models based 
on the Holstein Hamiltonian, the effects of the phonon 
spectrum renormalization on the polaron scattering have 
been studied both for the optical [18] and the acoustical 
[19] branches. Phonon anharmonicities have been also 
proposed as a fundamental feature favouring high T c su- 
perconductivity in the context of (bi)polaronic theories 
[20-22]. 

When the electron energy scale, set by the value of 
the tight binding overlap integral, is larger of (or com- 
parable to) the characteristic phonon energy the lattice 
deformation does not follow instantaneously the electron 
motion [23] and the e-ph interaction becomes time de- 
pendent. The path integral formalism [24] is suitable 
to account for time dependent interactions as a retarded 
potential naturally emerges in the exact integral action 
[25]. Moreover, being valid for any e-ph coupling value, 
the path integral method allows one to derive the par- 
tition function and the related temperature derivatives 
without those limitations which affect the perturbative 
studies. 

Our previous investigations have focussed on the equi- 
librium thermodynamics of the SSH model by assuming 
a bath of harmonic phonons for the electron particle path 
[26] while anharmonic corrections to the phonon free en- 
ergy had been only partially considered [27] in the case 
of a time averaged electron particle path. However, be- 
ing the e-ph SSH interaction intrinsecally non local on 
the time scale, significant phonon anharmonicities may 
be found through a full computation of the cumulant 
expansion depending on the perturbing source current 
of the Hamiltonian model. With the present paper we 
propose a comparative analysis of the equilibrium ther- 
modynamics for the phonon system in the Holstein and 
SSH models pointing out the appearance of macroscopic 
anharmonic features driven by the strength of the e-ph 
coupling. In Section II we outline the Hamiltonian mod- 
els and the cumulant expansion for the phonon partition 
function. The results are presented in Section III, while 
Section IV contains some final remarks. 



II. THE HAMILTONIAN MODELS 

In one dimension the e-ph Holstein Hamiltonian (H H ) 
and the e-ph SSH Hamiltonian (H SSH ) read respectively: 

H " = J 2 £(/r/r+l + fl + ,fr) + s£ + W) 

r r 

jjSSH _ J r r+1 (/t/ r+1 + / r+1 / r ) 
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Jr,r+1 = — ^ [J + a ( u r - Ur+l)] 

(0.1) 

where, 5 is the Holstein e-p/i coupling in energy units, 
a is the SSH e-p/i coupling in units energy x length -1 , 
J is the nearest neighbors hopping integral for an undis- 
torted chain, /t and / r create and destroy electrons on 
the r— lattice site, u r is the displacement of the r — th 
atomic group along the molecular axis, b I and b r are 
the real space creation and annihilation phonon opera- 
tors. The non-interacting phonon Hamiltonian H p h is 
assumed to be given, in both models, by a set of in- 
dependent oscillators. Taking u> as the characteristic 
phonon energy, the displacement field can be written as 
u r = (2Mw)~ 1/2 (&t + b r ) and H H transforms into 

r r 

g = gV2Mu (0.2) 

with g dimensionally equivalent to the coupling a of 
the SSH model. Choosing the atomic mass M of or- 
der 10 4 times the electron mass, we get g ~ 1.1456 x 
g\f2u) \meVA~ 1 ] where g is given in units of u> and the 
latter is given in meV. This transformation allows us 
to compare numerically the e-ph effective coupling in the 
two Hamiltonian models. 

Let's apply to the Hamiltonians in eq.(l) space-time 
mapping techniques previously used i.e. in the path in- 
tegral formulation of the theory of dilute magnetic alloys 
[28] and in the path integral study of A15 compounds 
with strong electron-phonon coupling [29]. 

Thus, we introduce x(t) and y(r') as the electron co- 
ordinates at the r and r + 1 lattice sites respectively, 
while the spatial e-ph correlations contained in (1) are 
mapped onto the time axis by changing: u r — ► u{t) and 
u r+1 -> u(r'). Then, H H in (2) and H SSH in (1) trans- 
form to 

H h (t,t')= ~(p(x(r))f(y(T')) + p(y(T'))f(x(r))) 

+ gu{r)p{x{r))f{x{r)) 
H ssh (t,t') = J T , T ,(p(x(T))f(y(T>)) + f\y(r'))f(x(r))) 

J T , T > = ~\[J + a{u{r) - u{t'))} (0.3) 

r and t' are continuous variables (<G [0,0]) in the Mat- 
subara Green's functions formalism with being the in- 
verse temperature hence, the electron hops are not con- 
strained to first neighbors sites along the chain. Accord- 
ingly, the Hamiltonians in (3) are more general than the 
real space Hamiltonians in (1). Note that, while the Hol- 
stein model assumes a local e-ph interaction along the t— 
axis, in the SSH picture the electron hopping is accom- 
panied by a time dependent displacement of the atomic 
coordinate. 
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The ground state of the SSH Hamiltonian is twofold 
degenerate and, in real space, a soliton connects the two 
phases with different senses of dimerization. As a local- 
ized electronic state is associated to the soliton the SSH 
model describes in principle both electron hopping be- 
tween solitons and thermally activated hopping to band 
states. By mapping onto the t— axis we introduce time 
dependent electron hops while maintaining the funda- 
mental features of the SSH Hamiltonian. Eqs.(3) display 
the semiclassical nature of our method in which quantum 
mechanical degrees of freedom interact with the classical 
variables u(t). After setting r' = 0, u(0) = y(0) = 0, we 
take the thermal averages for the electron operators over 
the ground states of the respective Hamiltonians thus ob- 
taining the average energies per lattice site due to elec- 
tron hopping plus e-ph coupling: 



< H h (t) > _ J 



N 

< H ssh (t) > 



-G±[x(t),t] + gu(r) 



N 



-[- + au(r))G ± [x(r),r} 
G ± [x(t),t] = G[-x(t),-t]+G[x(t),t] 

= — / dk cos{k x(t)] cosh (ekT)np(ek) 
77 Jo 

(0.4) 

where N is the number of lattice sites, G[x(t),t] is 
the electron propagator, a is the lattice constant, Ek = 
— t; cos(fca) and np{£k) is the Fermi function. 

In general, the phonon partition function perturbed 
by a source current j(r) can be expanded in anharmonic 
series as: 



Z ph[j(T)} 
N 

z h =U 



zji+j2(-iy<c i > j{T) ) 
\ i=i * 



2sinh(wi/3/2) 



(0.5) 



where the cumulant terms < C l >j< T ) ar e expectation 
values of powers of correlation functions of the perturb- 
ing current. The averages are meant over the ensem- 
ble of the N harmonic oscillators (having energies u>i) 
whose partition function is Z^. As the oscillators are de- 
coupled we study the anharmonicity induced on a single 
oscillator with energy u> by the source currents peculiar 
of the Hamiltonians in (3). Being, j H (t) = gu(r) and 



,ssh {t) 



< C k > jH = 



■au(r)G ± [x(t), t] we get 
* Du{t) 



k\ J ~ 



xexp 



< C k > jS SH (x) = 



f 

Jo 



k\ 



/ cLtu(t) 
Jo 



dT—{u 2 {T) + 



2 2 
UJ U 



(r)) 



j> Du{t) 






Jo 



dTu(T)G ± [x{T),T] 
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xexp 



dT^-(u 2 (T)+u; 2 u 2 (r)) 



As a direct consequence of the time retardation in the 
SSH e-ph interactions, the SSH cumulant series turns out 
to depend on the electron path coordinates. 

In order to calculate eqs.(6), we expand the periodic 
oscillator path (u(t + (3) = u(t)) in Np Fourier compo- 
nents: 



u{t) = u a + ^ 2($iu n cos(w„r) — 3u„ sin(o> rl r)^ 



2im/f3 



(0.7) 



and choose a measure of integration which normalizes 
the kinetic term in the oscillator field action: 



<j> Du{t) 



V2 [°° 

, s2N F + l / dU ° 
N p 



\{(2nnf 



Du{r)exp 



r M f (: 



dStu n 



dru (r) 



d$su n 



(0.8) 



being Am — 

As Jq dTu(r) = /3uo, using eqs.(7) and (8), the Hol- 
stein cumulant series can be worked out analytically as 
a function of the cutoff Np. The first of eqs.(6) gets: 



< C fc > 



k ^N B 



z h -*w\ M ) k (k-\)\\ g 



(2n7r) 2 



7r fe / 2 (w/3) 



fe+i 



\ (2n7r) 2 + (ujl3) 2 
(0.9) 



where only even k terms survive. As the cumulants 
should be stable against the number of Fourier compo- 
nents in the oscillator path expansion we set the mini- 
mum Np through the condition 2Npir ^> ujj3. Since the 
Holstein source current does not bear any dependence 
on the electron path coordinates, the Holstein expansion 
in eq.(9) coincides with the cumulant series for a time 
averaged electron particle path [27]. 

To compute the second of eqs.(6), we expand also the 
electron path under the closure condition x(t + /?) = 
x(t): 

M F 

x(t) = x + ^ 2^Rx m cos(w m r) - 3x m sin(w m r)^ 



u) m = 2irm/(3 



(0.10) 



with, in general, Mp ^ Np. Thus, the general cumu- 
lant term is obtained by averaging the second of eqs.(6) 
over a set of electron particle paths: 



(0.6) 
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S S H = SSH (x) (0.11) 

Two Fourier components (Mp = 2) suffice to get nu- 
merical convergence in eq.(ll). We point out that, the 
computation of the anharmonic expansion in the SSH 
model requires for any choice of electron path coefficients, 
double summations over Brillouin zone (see eq.(4)) and 
inverse temperature (r) values (see the second of eqs.(6)). 



III. PHONON THERMODYNAMICS 

The phonon thermodynamics can be derived from 
eq.(5) via computation of the phonon free energy 
FW(T)= -(5-Hn[Z ph \j{r)]]. " 

To proceed one needs a criterion to find the tempera- 
ture dependent cutoff k* in the cumulant series. We feel 
that, in the low T limit, the third law of thermodynamics 
may offer the suitable constraint to determine k* . Then, 
set the input parameters a (g in the Holstein model), J 
and w, the program searches for the cumulant order such 
that phonon entropy Ey -* (T) and heat capacity Cy •* (T) 
vanish in the zero temperature limit. In a constant vol- 
ume transformation, we compute 

E { y ] (T) = - (T + 2 A) - 2F™ (T + A)] / A 

(T) = - (T + 2 A) - 2F^ (T + A) + (T) 

x (k + z&) (0 - 12) 

A being the incremental step and k* is determined, at 
finite T, as the minimum value for which the anharmonic 
cumulant series converges with an accuracy of 10~ 4 . The 
results are displayed in Figures 1 for the case of a low 
energy oscillator and a narrow band hopping integral, 
J = O.leV. Three values of a are assumed to emphasize 
the dependence of the SSH model thermodynamics on 
the e-ph coupling. The g = 3 value for the Holstein 
model calculations yields a g coincident with the largest 
a = 21.74meV r l" 1 . 

The Holstein cumulant expansion causes a slight and 
slope preserving lowering of the phonon free energy with 
respect to the harmonic plot. As a consequence the Hol- 
stein entropy and heat capacity, computed via eqs.(12), 
are superimposed to the harmonic quantities. However, 
looking at Fig. 1(d) (case g = 3), we see that a large num- 
ber of cumulant terms is required to get convergence in 
the Holstein series as the cutoff k* markedly grows below 
SOA', attaining the value k* = 40 at T = 1 A. Thus, e-ph 
anharmonicities strongly renormalize the Holstein parti- 
tion function, essentially stabilize the system, but there 
is no trace of anharmonicity in the phonon free energy 
temperature derivatives [30]. 

Much different is the trend of the anharmonic correc- 
tions in the SSH model. The main features can be sum- 
marized as follows: i) the phonon free energy (Fig. 1(a)) 
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and its derivatives (Figs. 1(b), (c)) largely deviate from 
the harmonic plots and strongly depend on a; ii) there 
is a threshold a value above which the heat capacity dis- 
plays a peak. By enhancing a (Fig. 1(c)) the height of the 
peak grows and the bulk of the anharmonic effects on the 
heat capacity is shifted towards lower T; iii) at high T 
the anharmonic corrections renormalize downwards the 
phonon free energy but their effect on the heat capac- 
ity tends to decrease signalling that e-ph nonlinearities 
are rather to be seen in the intermediate temperature 
range. Note that the size of the anharmonic enhance- 
ment is about a factor 4 larger than the value of the 
harmonic oscillator heat capacity at T ~ 170 K; iv) the 
cutoff k* in the SSH model (Fig. 1(d)) attains the value 
10 at room temperature while the cumulant series con- 
verges with a low number of terms (k* = 2) at low T. 
At intermediate T £ [90, 270] , the SSH cutoff overlaps 
the Holstein cutoff at the value k* — 8. 

The location of the heat capacity peak along the 
temperature axis strongly varies with the oscillator en- 
ergy. In Fig. 2 we select two energy values and report 
the crossover temperatures for a set of e-ph coupling 
strengths. By increasing ui at fixed a the peak shifts 
towards high T and it shows up above a minimum a 
whose value increases at larger u. Thus, higher energy 
oscillators are more stable against e-ph induced anhar- 
monicity since larger couplings are required to introduce 
signatures of non linear behavior in the phonon heat ca- 
pacity. 

So far we have focussed on the thermodynamics of the 
phonon subsystem. The question arises whether and to 
which extent the anharmonic phonon heat capacity af- 
fects the total heat capacity of the SSH model. As dis- 
cussed in ref. [26] the total heat capacity is given by the 
phonon contribution plus a source heat capacity which 
includes both the electronic contribution (related to the 
electron hopping integral) and the contribution due to 
the source action (the latter being oc a 2 , see eqs.(6) in 
ref. [26]). A bath of harmonic phonons has been as- 
sumed in our previous work. Now we compare the an- 
harmonic phonon heat capacity (C v nh ) and the source 
heat capacity, here termed C v ~ p to emphasize the de- 
pendence both on the electronic and on the e-ph coupling 
terms. The plot displaying the largest peak in Fig. 1(c) 
is reported on in Fig. 3 (a) while C' v p is computed ac- 
cording to the method described in in ref. [26] after set- 
ting a = 21.74meVA _1 . Also the total heat capacity 
(C^ ot = C$ nh + C e y p ) is shown in Fig.3(a). At low tem- 
peratures, Cy~ p yields the largest effect mainly due to 
the electronic hopping while at high T, C^ p prevails 
as the source action becomes dominant. In the inter- 
mediate range (T £ [90,210]) the anharmonic phonons 
provide the highest contribution although their charac- 
teristic peak is substantially smeared in the total heat ca- 
pacity by the source term background. Fig. 3(b) compares 
Cy nh and C v ~ p for two increasing values of a: while the 
anharmonic peak shifts downwards (along the T axis) by 
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enhancing a, Cy nh remains larger than C^T P in a temper- 
ature range which progressively shrinks due to the strong 
dependence of the source action on the strength of the 
e-ph coupling. Finally, it has to be pointed out that the 
low temperature upturn displayed in the total heat ca- 
pacity over T ratio [26] is not affected by the inclusion 
of phonon anharmonic effects which tend to become neg- 
ligible at low temperatures. 

IV. FINAL REMARKS 

The Holstein and the Su-Schricffcr-Heeger Hamilto- 
nian present a fundamental difference with respect to the 
nature of the e-ph coupling. Mapping the real space in- 
teractions onto the time scale we have shown that the 
Holstein e-ph interactions are local in time while, in the 
SSH model, the electronic hopping induces a time depen- 
dent lattice displacement thus leading to a non local e-ph 
coupling. As a consequence the source current peculiar 
of the SSH model depends both on time and on the elec- 
tron path coordinates. We have applied the path integral 
method to analyse the perturbing effects of the source 
currents (in both models) on the phonon subsystems by 
expanding the phonon partition function in cumulant se- 
ries. In the low temperature limit, the constraint imposed 
by the third law of thermodynamics permits to deter- 
mine the cutoff k* in the anharmonic expansion: while 
the SSH series converges rapidly, an high number of cu- 
mulants is required in the Holstein model to stabilize the 
low T equilibrium properties. At intermediate and high 
T, k*(T) is comparable in both models. However, look- 
ing at the behavior of the thermodynamic properties we 
find striking differences between the two models. The 
Holstein phonon heat capacity does not show any trace 
of anharmonicity induced by the e-ph coupling while the 
SSH phonon heat capacity exhibits a peak whose location 
along the T axis strongly depends on the strength of the 
coupling a. By increasing a at a fixed T, k* grows and 
the phonon partition function becomes larger. As a result 
the point of most rapid decrease in the SSH phonon free 
energy versus temperature shifts downwards and the as- 
sociated heat capacity peak is found at lower T. At high 
temperatures the free energy decreases regularly and the 
anharmonic effects on the heat capacity are reduced. In- 
fact, at high T, the time (r) retardation between electron 
hopping and atomic displacement is shortened as r en- 
counters the upper limit of a small inverse temperature. 

Finally, we have computed the general path integral in 
the SSH model for an electron path moving in a bath 
of anharmonic oscillators. We have derived the total 
heat capacity and pointed out the contribution of the 
anharmonic phonon with respect to the purely electronic 
and e-ph terms. The anharmonic effects are relevant in 
the intermediate temperature range whereas the peak is 
smeared by the electronic hopping term (on the low tem- 
perature side) and by the source action heat capacity at 
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high temperatures. 



FIG. 1. Phonon (a) Free Energy, (b) Entropy and (c) 
Heat Capacity versus temperature in the Su-Schrieffer-Heeger 
model with e-ph anharmonic contributions due to three val- 
ues of the coupling a in units meV A~ l . The Holstein model 
results are computed for the coupling g (in units of ui) cor- 
responding to the largest of the three a values. Also the 
harmonic phonon functions are plotted for comparison, (d) 
Number of anharmonic terms in the cumulant series yielding 
convergent thermodynamical quantities versus temperature 
both in the Su-Schrieffer-Heeger model (a case) and in the 
Holstein model (g case). A low energy oscillator is assumed. 
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FIG. 2. Temperature locations of the peak in the 
Su-Schrieffer-Heeger phonon heat capacity versus e-ph cou- 
pling in units meVA^ 1 . 



FIG. 3. (a) Total Heat Capacity versus temperature in the 
Su-Schrieffer-Heeger model. The contributions due to anhar- 
monic phonons (C v nh ) and electrons plus electron-phonon in- 
teractions (Cy~ p ) are plotted separately. The largest a of 
Fig.l is assumed, (b) Cy nh and C|T P for two values of a in 
units meVA~. ui — 20meV. 
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